Search for: All records

Documenting the processes and timescales of magma formation and diversification and defining the locations, shapes, volumes, and phase states of magma storage and transport zones rely on data produced by novel analytical techniques and state-of-the-art experimental methods. Computational modeling effectively links these critical tools. The Magma Chamber Simulator (MCS) is an internally consistent thermodynamic open system model that uses experimental constraints from rhyolite MELTS (Gualda et al. 2012, Ghiorso & Gualda, 2015) to compute paths of open system magmas that evolve via processes including crystallization, magma mixing, cumulate/mush entrainment, and host-rock assimilation. MCS results yield elemental, isotope, mass, and thermal characteristics of melt ± crystals ± volatiles in "resident" magma, crustal wallrock (melt and solids), recharge magma, and entrained material. To model the petrochemical evolution of igneous rocks related by open-system processes, one typically runs 200+ models that vary initial compositions, pressures, and temperatures of magma, wallrock, etc. Comparison of model results with whole rock, mineral, and melt inclusion chemical data and other constraints (e.g., thermobarometry) yield interpretations about igneous processes at a range of scales—from how crust forms and evolves to processes responsible for in situ geochemical records of crystals—and allows assessment of epistemic and aleatoric uncertainties. Two examples of computational studies will illustrate MCS's utility and flexibility. (1) Modeling of historical basalts at Mt. Etna (Italy) provides evidence for variable degrees of melting of metasomatized mantle, followed by magma recharge and assimilation of partial melts of carbonate-flysch crust. (2) MCS models reproduce whole rock and mineral data of plagioclase-rich basalts at Steens Mountain (USA) through entrainment of gabbroic mush that likely formed in early stages of Columbia River Basalt magmatism. To enhance understanding of trans-lithospheric magma systems, future work on MCS will prioritize (i) building a post-processing environment that utilizes select statistical methods to inform "best-fit" models and to quantitatively assess uncertainty, and (ii) increasing modeling efficiency by adding automated modeling capabilities. 

Abstract The nearly continuous volcanic eruption record at Mt. Etna dating back ~700 years provides an excellent opportunity to investigate the geochemical evolution of a highly active volcano. Of particular interest is elucidating the cause of selective enrichment in alkali elements (K and Rb) and 87Sr/86Sr observed in various episodes of past activity. More recently, this alkali enrichment trend started to manifest in the 17th century and accelerated after 1971, and was accompanied by an increase in the volume, frequency, and explosivity of eruptions. Explanations for this signature include recharge of alkali-enriched magmas and/or crustal contamination from the subvolcanic basement. This study quantitatively examines the role of crustal contamination in post-1971 Etnean magma compositions via hundreds of open-system phase equilibria and trace element calculations based upon whole-rock major oxides, trace elements, 87Sr/86Sr ratios, and mineral compositional data. Available pre-1971 petrochemical data are satisfactorily reproduced by fractional crystallization of a high whole-rock MgO (12–17 wt.%), Ni (135–285 ppm), and Cr (920–1330 ppm) parental magma composition that is documented in Etna's ~4-ka fall-stratified deposit. Observed post-1971 whole-rock and glass trends and phase equilibria are reproduced via modeled assimilation of a skarn and flysch mixture, lithologies that represent the uppermost 10 to 15 km of sedimentary rocks beneath Etna. Notably, models show that K2O (wt.%) and Rb (ppm) behave incompatibly during partial melting of skarn/flysch. Additionally, the observed elevation of 87Sr/86Sr in post-1971 samples is consistent with the addition of radiogenic Sr from wallrock partial melts. In best-fit models, which yield observed post-1971 K2O, Rb, and 87Sr/86Sr trends, ~17% anatectic melt is assimilated and there may be a subordinate stoped wallrock component of ≤2% (percentage is relative to the starting mass of pristine magma). Previous work has shown that metasomatized spinel lherzolite and garnet pyroxenite can be melted in different proportions to reproduce long- and short-term changes observed in Etna’s geochemical products. We propose that the alkali enrichment signature observed after 1971 can be fully explained through the combination of mantle heterogeneity and crustal contamination. In particular, up to ~20% crustal input coupled with mantle heterogeneity of primitive melts explains the geochemical signals quite well. The influence of crustal contamination on post-1971 lavas is, in part, the result of frequent recharge of magmas that thermally primed the middle to upper crust and enhanced its partial melting.